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ABSTRACT 

Aims. Solar convection in strong plage, having a magnetic field that is vertical and strong over extended regions, but much weaker 
than in umbrae of large sunspots, has so far not been well studied. This is mostly due to lack of spectropolarimetric data at adequate 
spatial resolution. The combination of a large solar telescope, such as the Swedish 1-m Solar Telescope, adaptive optics, powerful 
image reconstruction techniques and a high-fidelity imaging spectropolarimeter is capable of producing such data. In this work, we 

I ^ study and quantify properties of strong-field small-scale convection and compare observed properties with those predicted by numer- 

ical simulations. 

r^ Methods. We analyze spectropolarimetric 630.25 nm data from a unipolar ephemeral region near sun center We use line-of-sight 

Cn velocities and magnetic field measurements obtained with Milne-Eddington inversion techniques along with measured continuum 

intensities and Stokes V amplitude asymmetry at a spatial resolution of 0.'15 to establish statistical relations between the measured 

■j' quantities. We also study these properties for diff'erent types of distinct magnetic features, such as micropores, bright points, ribbons, 

Ph flowers and strings. 

^0 Results. We present the first direct observations of a small-scale granular magneto-convection pattern within extended regions of 

strong (more than 600 G average) magnetic field. Along the boundaries of the flux concentrations we see mostly downflows and 
asymmetric Stokes V profiles, consistent with synthetic line profiles calculated from MHD simulations. We note the frequent occur- 
rence of bright downflows along these boundaries. In the interior of the flux concentrations, we observe an up/down flow pattern 
Q that we identify as small-scale magnetoconvection, appearing similar to that of field-free granulation but with scales 4 times smaller. 

?-H Measured RMS velocities are 70% of those of nearby field-free granulation, even though the average radiative flux is not reduced. 

^ The interiors of these flux concentrations are dominated by upflows. 

I I Key words. Convection - Sun: faculae, plages - Sun: granulation - Sun: Magnetic fields - Polarization 

^ 1. Introduction Solar Telescope (SST; Scharmer et al. 2003a) combined with 



Oh 



l> 



3D MHD simulations (Keller et al. 2004; Carlsson et al. 2004; 



^O The interaction of convection with an existing initially weak 

^ magnetic field is the primary step in the production of the strong In contrast to bright small-scale magnetic features, the 

|> small-scale flux concentration features observed in the solar pho- largest magnetic flux concentrations, sunspots, inhibit most of 

O tosphere. When the average field strength is small enough, the the convective energy flux from below and lead to intermittent 

O horizontal flow component of convection sweeps the magnetic i^^arly field-free convection in the form of narrow plumes, ob- 

T^H field to the intergranular lanes, leading to expulsion of field from served as umbral dots, when the field is strong and nearly verti- 

^ the center of granules and flux concentrations at the intergranular ^^1 (Schiissler & Vogler 2006). 

• ^-H lanes. These flux concentrations reduce convective heat transfer, In the present paper, we discuss observations of magnetic 

/\^ leading to super adiabatic cooling and collapse into a strong- structures with scales corresponding to bright points and larger 

^ field state (Parker 1978; Spruit 1979), referred to traditionally butwithfieldstrengthslowerthanthoseof sunspots. These small 

as a flux tube or, when more appropriate, as a flux sheet. Once scale structures include abnormal granulation (Dunn & Zirker 

formed, a flux tube remains in a strong-field state for some time, 1973; Title et al. 1992), micropores and other distinct small-scale 

even though it evolves continuously as a result of surrounding magnetic features that have been dubbed ribbons, flowers and 

granular convection. If the filling factor of these flux tubes is strings (Berger et al. 2004). The impact of the magnetic field 

sufl&ciently small, the result is a granular pattern that is rela- on the convection within and adjacent to these magnetic field 

lively unaflTected by the flux tubes but with adjacent downflows concentrations has so far not been well studied, mostly due to 

driven by enhanced radiative cooling through the optically thin lack of adequate data at sufl&cient spatial resolution. By means of 

flux tubes. The appearance of flux tubes, in particular when ob- spectropolarimetric data and inversion techniques, we describe 

served away from disk center or in spectral lines, are as "bright convective signatures separately for abnormal granulation and 

points" or faculae when the flux concentrations are larger. The the various types of small-scale magnetic structures described in 

formation and evolution of such isolated flux concentrations are Berger et al. (2004) and compare these observations to existing 

now well understood from observations using the Swedish 1-m 3D MHD simulations. 
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2. Observations and data reduction 

2.1. Instrumentation 

The observational data were recorded with the SST and the 
CRisp Imaging SPectropolarimeter (CRISP; Scharmer et al. 
2008). CRISP is a telecentric dual Fabry-Perot etalon sys- 
tem used in combination with two nematic liquid crystal vari- 
able retarders (LCVRs) and a polarizing beam splitter. CRISP 
uses a unique combination of a high-resolution, high-reflectivity 
etalon with a low-resolution low -reflectivity etalon. This pro- 
duces a high-fidelity transmission profile with high transmission 
and small variations of its shape over the field-of-view (FOV) 
(Scharmer 2006). The polarizing beam splitter separates the 
beam from CRISP into two orthogonally polarized states, allow- 
ing strong suppression of seeing-induced cross-talk from Stokes 
ItoQ,U and V. The two linearly polarized images are recorded 
with separate CCD's and cover about 70''x70'' at an image scale 
of Off 01 1 /pixel (the image scale was changed to O'.'059/pixel in 
2009). A third CCD is used to record broad-band images through 
the prefilter of CRISP, providing needed support for reconstruc- 
tion and alignment of the narrowband images obtained at dif- 
ferent wavelengths and polarization states (see below). Images 
recorded through the three CCD's are exposed simultaneously 
by means of a rotating chopper. To improve image quality and 
image restoration, short exposure times are used. The frame rate 
was set by the chopper at 36 Hz, corresponding to an exposure 
time of about 17 ms and a "dark" (CCD readout) time of 1 1 ms. 
To compensate for the relatively poor signal-to-noise of each 
short exposure frame, a large number of frames are needed at 
each wavelength and polarization state. 



2.2. Observations 

We observed a unipolar (within the FOV) ephemeral region, lo- 
cated close to disk center (heliocentric coordinates S9 Ell) on 
26 May 2008. CRISP was used to record images at 9 wave- 
lengths between -19.2 and +19.2 pm in steps of 4.8 pm from the 
line center of the Fe 1 630.25 nm line and also at an adjacent con- 
tinuum wavelength. At each wavelength position images were 
obtained with 4 diff'erent LCVR tunings to allow measurements 
of the complete Stokes vector. For each wavelength and LCVR 
state, 16 images were recorded per camera. Each sequence pro- 
cessed thus consists of up to 640 images per CCD (1920 images 
in total) but in practice missing or corrupted images reduced the 
number of useful frames per camera to about 630. Each data set 
corresponds to a wall clock time of approximately 21s, includ- 
ing time for tuning CRISP. A time series of about 40 minutes was 
recorded in reasonably good seeing conditions during which the 
SST adaptive optics system (Scharmer et al. 2003b) indicated a 
lock rate between 50% and 90% most of the time. In this paper, 
we discuss only one of the best snapshots recorded. 



2.3. Data reduction 

The images from each 21 sec scan were processed as a sin- 
gle MOMFBD data set using the image reconstruction method 
developed by Lofdahl (2002) and implemented by van Noort 
et al. (2005). In addition to compensating residual seeing par- 
tially compensated by the AO system, this processing produces 
restored images at diff'erent wavelengths and polarization states 
that are accurately co-aligned with respect to each other. The 
processed images were demodulated and corrected for telescope 
polarization, using the telescope polarization model of Selbing 



(2005). This model typically gives residual cross-talks, or errors 
in the Miiller matrix elements, that are on the order of 1%, the 
dominant eff'ect of which is cross-talk from Stokes / to Q, U, 
and V. This cross-talk was compensated for with the aid of the 
Stokes images recorded in the continuum. There is most likely 
also cross-talk from V, contaminating the Q and U profiles (and 
vice- versa). However, due to the proximity of this region to sun 
center, the Q and U profiles are in any case very weak and nei- 
ther provide meaningful information about field inclinations nor 
field azimuths. We therefore did not attempt to compensate for 
cross-talk from V to Q and U and in the following discuss only 
the Stokes V profiles and the associated LOS component of the 
magnetic field. 

To visualize (by enhanced contrast) the various small-scale 
magnetic structures discussed in Sect. 3.3, we also calculated 
maps of the 630.25 nm line minimum intensity. This was done 
by fitting the center part of the Stokes / line profiles to a second- 
order polynomial and from the coefficients of these fits calculat- 
ing the minimum intensity at each pixel. 



2.4. Inversions 

Inversions of the data were carried out using the Milne- 
Eddington inversion code Helix (Lagg et al. 2004). The inver- 
sions involve calculation of the synthetic Stokes profiles based 
on Unno-Rachkovsky solutions of the polarized radiative trans- 
fer equations with the fitting and optimization controlled by the 
genetic algorithm Pikaia (Charbonneau 1995). To account for 
the limited spectral resolution of CRISP, about 6.4 pm at 630 nm, 
the synthetic profiles were convolved with the theoretical CRISP 
spectral transmission profile at each iteration. The following free 
parameters were solved for with Helix: The field strength B, the 
azimuth angle, the inclination angle y, the line-of-sight (LOS) 
velocity vlos, the Doppler width and the gradient of the source 
function. The fixed parameters were the line- strength (ratio of 
line to continuum opacity), set to 16, and the damping param- 
eter, set to 1. This choice was made based on numerous tests 
with Helix using the data and inversions discussed by Scharmer 
et al. (2008). Fits were made with diff'erent weights for the 4 
components of the Stokes vector. Following the suggestion of 
A. Lagg (private communication), we used weights for / and V 
that were twice as large as for Q and U. Extensive tests with 
Helix and the CRISP data discussed by Scharmer et al. (2008) 
failed to produce stable estimates of the magnetic filling factor 
/, but demonstrated that / multiplied by the field strength was 
consistently close to what was obtained by setting / to unity. 
Also for the present inversions, we therefore set / to unity, al- 
though already the limited spatial resolution implies the exis- 
tence of spatial stray light for unresolved magnetic structures. In 
addition, a diff'use component of straylight must exist to explain 
the measured granulation contrast (see Sect. 2.6). The measured 
LOS component of the magnetic field ^los. discussed in the fol- 
lowing sections of the paper, therefore corresponds to a magnetic 
flux density or a (poorly defined) spatial average. This measured 
average most likely is an underestimate in magnitude near the 
centers of magnetic flux concentrations and probably an overes- 
timate at or outside their magnetic boundaries because of stray- 
light. Due to the proximity of the observed region to sun center 
and the absence of systematic difl'erences between the disk cen- 
ter and limb side of the observed quantities discussed, we did 
not transform the measured magnetic field vectors to the solar 
frame. 
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Fig. 1: Full Field of view showing the unipolar ephemeral region with small pores with the 2r'x27''region of interest (ROI) marked 
as a red or black rectangle. Tick marks are in units of arcsec. The upper row shows the restored CRISP continuum image.The other 
three panels show quantities returned by the inversion code: the LOS velocity (lower-left) and LOS magnetic field to enhance the 
strong (upper-right) and weak (lower-left) field. The LOS velocities are shown enhanced by adding an off'set to the velocities within 
the 200 G contour. The grey level bars on the right hand side of the figure indicate the signed LOS magnetic field in units of Gauss. 



2.5. Velocity calibration 

To compensate for the wavelength shifts due to cavity errors 
dominated by the high-resolution etalon, we used the flat-field 
images, obtained 1 hour after recording the science data. At each 
pixel, we made a fit of a second-order polynomial to the cen- 
ter part of the 630.2 nm line profile and measured the wave- 
length shifts. These wavelength shifts (converted to equivalent 
Doppler shifts) were subtracted from the LOS velocities returned 
by Helix. The zero point of the corrected LOS velocities was 
fixed by assuming the dark parts of the pores to be at rest and us- 
ing a mask based on the restored continuum image to identify the 
pores. Adopting this zero point and calculating the average LOS 



velocity outside the center 400x400 pixels, we obtained a con- 
vective blueshift of -270 ms"^ The obtained blueshift is close 
to the convective blueshift of -210 m s"^ obtained by Scharmer 
et al. (2008), based on a similar set of CRISP data also recorded 
close to sun center and also using the Helix inversion code. We 
also compared to convective blueshifts obtained from 3D con- 
vection simulations (de la Cruz Rodriguez, private communica- 
tion). An average 630.2 nm line profile calculated from these 
simulations produces a bisector that midways between the con- 
tinuum and the deepest part is shifted by about 300 m s"\ which 
is close to the value we obtained. This suggests that our velocity 
calibration is accurate to within about 100 m s"^ but the obtained 
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convective blueshift appears on the high side and we cannot ex- 
clude an error of 200 m s"^ A more accurate assessment would 
require inversions with Helix using synthetic line profiles calcu- 
lated from the convection simulations and degraded to the spec- 
tral resolution and sampling of the CRISP data. 



2.6. Data quality and fits 

The polarimetric noise was estimated from a 500 x 500 pixel 
area in the Stokes Q, U and V images recorded at the contin- 
uum wavelength and corresponds to 1.2 x 10"^, 1.4 x 10"^, and 
1 .Ox 10"^ for Q, U and V respectively. The continuum Stokes im- 
ages appear completely clean from granulation features or other 
artifacts, suggesting very low levels of seeing induced cross-talk. 
We also measured the RMS noise in ^los over several 20x20 
pixel boxes that appear free of polarization features and obtained 
noise levels in the range 2.2-3.2 G, suggesting "safe" detection 
levels of about 10-20 G. 

Figure 1 shows a map of ^los scaled from -1000 G to 
+370 G (upper right panel) and from -20 G to +20 G (lower 
panel). When scaled to show weak field, the center part of the 
FOV appears to have magnetic field nearly everywhere. This 
may be in part straylight from uncompensated high-order aber- 
rations in MOMFBD processing or other sources of straylight 
in the telescope and/or the Earth's atmosphere. The presence 
of straylight is clearly indicated by the measured RMS contrast 
from the CRISP continuum image. The high continuum RMS 
contrast measured, about 8.3% in the region without significant 
magnetic field, is still far from the value obtained from 3D con- 
vection simulations, 14%. This is not a spatial resolution issue: 
the main contributions to the granulation RMS contrast come 
from spatial frequencies that are well below the SST diff'rac- 
tion limit. We also note that the strongest LOS component of the 
magnetic field shown in Fig. 1 is only about 1.2 kG which ap- 
pears low. Including the transverse magnetic field components, 
the highest field strength is about 1.3 kG. Inversions based on 
Hinode spectropolarimetric data from a similar region with small 
pores (Morinaga et al. 2008) gave peak mean field strengths 
(magnetic filling factor multiplied by the field strength) of 1.5- 
1.6 kG. 

Most observed V profiles, some of which are discussed in 
Sect. 3.3, are rather symmetric, justifying the use of a Milne- 
Eddington inversion technique to analyze the data. Only Stokes 

V profiles with very low amplitudes are poorly fitted. We also 
note a few cases where the shapes of the observed Stokes / and 

V profiles are well reproduced by the fits, but where the (very 
weak) observed Stokes V profiles are strongly redshifted with re- 
spect to the synthetic Stokes V profiles. This implies that the ob- 
tained LOS velocities primarily refer to the non-magnetic rather 
than the magnetic atmosphere or are contaminated with spatial 
straylight, where a diff'erence exists. The bias toward estimates 
of LOS velocities from the Stokes / rather than the V profiles is 
not surprising, given the equal weights on Stokes / and V used 
with the inversions and the very weak Stokes V signatures for 
these features. 

There is also a tendency for the synthetic Stokes / profiles 
to be somewhat redshifted with respect to the observed profiles. 
This is almost certainly due to the telluric blend in the red wing 
of the 630.2 nm line, which is not compensated for in the in- 
versions. However, there is a similar systematic influence of the 
telluric line also on the line profiles from the pores, used to cal- 
ibrate the LOS velocities, such that we expect the net eff'ect of 
the telluric blend on the measured LOS velocities to be small. 



3. Results 

3.1. Overview 

The restored CRISP continuum image is shown in the top panel 
of Fig. 1. The disk center (DC) direction is indicated with an 
arrow and the region of interest (ROI) discussed in Sect. 3.3 is 
indicated with a rectangle. Due to the proximity of this region to 
sun center, the transverse components of the magnetic field ob- 
tained with Helix are mostly weak and very noisy. In the follow- 
ing, we discuss the LOS component of the magnetic field, calcu- 
lated as ^Los = B cos y and scaled to show both the strong and 
weak fields in Fig. 1. In the MDI magnetograms this region is 
indicated as having positive polarity, but for greater clarity pos- 
itive polarities are here indicated as dark. The lower right panel 
demonstrates the high sensitivity of the data and the absence of 
spurious magnetic features. 

The observed region contains micropores in various stages 
of evolution and numerous flux concentrations, outlining a net- 
work of strong field located mostly in intergranular lanes. The 
region close to the pores contain more extended contiguous re- 
gions of strong magnetic field where the granulation pattern is 
small-scale and irregular, so-called abnormal granulation (Dunn 
& Zirker 1973). 

We note that the regions shown as dark in the lower-right 
panel of Fig. 1 contain no strong field of opposite polarity to 
that of this region and only a few cases with weak such field. 
Possibly, this is an eff'ect of polarized straylight from neighbor- 
ing strong field canceling the V signal from weak field of oppo- 
site polarity. 

3.2. Evidence for small-scale magneto-convection 

The lower left panel of Fig. 1 shows the LOS velocity obtained 
from the inversions, but with a mask used to enhance the re- 
gions where ^los is stronger than 200 G. The unsigned average 
of ^Los within this mask is 530 G, the average field strength is 
640 G. A mask with a threshold of 100 G (unsigned average of 
^Los 400 G, average field strength 510 G) or 300 G (unsigned 
average of ^los 620 G, average field strength 720 G) would out- 
line nearly the same regions but the 100 G mask would include 
also some of the larger granules embedded within the stronger 
field. We remark that the precise value of the threshold used 
to define the mask is of no particular physical significance for 
an extended region. Because of the finite spatial resolution and 
the critical sampling, any discontinuous change in the field will 
be observed as a continuous transition. The chosen threshold 
simply defines the boundary between the nearly field-free and 
strong-field regions as that for which I^losI is somewhat less 
than half the average of the adjacent strong field. 

There is a dramatic diff'erence in the velocity structures in- 
side and outside the mask. Within the 200 G mask (average field 
strength 640 G), we see a small-scale velocity pattern that ap- 
pears similar to the granular velocity field outside the mask, but 
with characteristic scales of about 0'.'3, or 4 times smaller than 
field- free granules. Also the continuum image shows distinctly 
small-scale, but more irregular, structures within the mask. The 
average LOS velocity within the mask is about 70 ms"^ The 
RMS velocity within the 200 G mask is 490 ms"^ compared 
to 700 ms"^ outside the mask. There appears to be some con- 
tributions to this RMS from large-scale velocity fields (see also 
lower-right panel in Fig. 6), possibly from 300-sec oscillations, 
implying that the RMS estimated for small-scale velocity field 
represents an upper limit. Using a 2"y2" boxcar average to re- 
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Fig. 2: The top row shows correlations between continuum intensity and LOS velocity for weak (I^losI <50 G), medium strong 
(200 G < I^LOsI < 800 G) and strong (800 G < I^losI) magnetic field and the correlation between the continuum intensity and LOS 
magnetic field. The bottom row shows correlations between LOS velocity and LOS magnetic field for dark (/c/</c) < 0.9), average 
(0.9 < Ic/ih) < 1.05) and bright (LOS < Ic/ih)) structures. Red lines show binned averages and in the lower row also ±1 standard 
deviation of the variations. 



duce the contributions from large-scale velocity fields reduces 
the RMS velocity within the 200 G mask to 450 m s"^ . We con- 
clude that the small-scale velocity pattern corresponds to up- 
ward/downward flows of relatively large amplitude that aver- 
ages approximately to zero. This RMS velocity decreases with 
increasing LOS field and is only 220 m s"\ if the 2''x2'' boxcar 
average is first subtracted, for the strongest fields measured. 

The RMS continuum intensity within the 200 G mask but ex- 
cluding the pores (defined by /c/</c) < 0.8, where </c) is the con- 
tinuum intensity averaged over the entire FOV) is 6.5% which is 
only 22% lower than that of the field-free granulation (8.3%). 
The top row of plots in Fig. 2 shows the correlation between 
continuum intensity and LOS velocity. The red curves indicate 
binned averages and ± one standard deviation of the variations. 
The left-most plot shows the correlation for nearly field-free con- 
vection (I^LOsI < 50 G), the second panel the same correlation 
with I^LOsI in the range 200-800 G and the third plot with I^losI 
stronger than 800 G. As shown in these plots, also the magnetic 
parts of the FOV show a (weak) correlation between the con- 
tinuum intensity and the LOS velocity in the same sense as for 
field- free convection. This correlation is stronger if the 2''x2'' 
boxcar average is used to reduce the eff'ects of large-scale flows 
and 300-sec oscillations, but such averaging was not used to pro- 
duce the plots shown. The relatively weak correlation found for 
the magnetic part of the FOV is obviously also related to the 
morphological difl'erences between the continuum intensity and 
LOS velocity maps; only the LOS velocity map shows a clear 
"granulation-like" pattern within the 200 G mask. These mor- 
phological difl'erences are very likely related to the formation 
height of the 630.25 nm line above the photosphere, the small 
horizontal scale of magnetic field variations inducing a similarly 
small vertical scale of such variations and the increasing domi- 
nance with height of the magnetic field, constraining the convec- 
tive flows above the photosphere. 

The small-scale LOS velocity pattern is barely visible at a 
resolution close to ~0'.'15. Quite obviously, the RMS velocities 
(as well as RMS intensity variations) measured within the mag- 



netic mask must be strongly underestimated, in particular since 
this is true already for the much larger field- free granulation pat- 
tern. The small-scale granular velocity pattern would be diffi- 
cult, if not impossible, to observe with a significantly smaller 
telescope. It is thus not surprising that Morinaga et al. (2008), 
using the 50-cm SOT on Hinode, did not report direct evidence 
for such small-scale convection, but interpreted their results as 
suppression of convection. 

Rimmele (2004) plotted the correlation between LOS ve- 
locities measured from two wavelengths in the wings of the 
630.25 nm and the line wing intensity for selected bright points 
(Rimmele 2004, Fig. 7). This correlation plot appears quite simi- 
lar to the first and third panel in the upper row of Fig. 2. Rimmele 
also plotted the correlation between ^los. estimated from single- 
wavelength magnetograms in the blue wing of the 630.25 nm 
line and the line wing intensity as well as the correlation be- 
tween vios and ^los (Rimmele 2004, Fig. 6). Also these plots 
appear similar to the corresponding plots in Fig. 2. 

We believe that the measured properties of the small-scale 
velocity field and the correlations shown in Fig. 2 justify the 
identification of the velocity pattern inside the mask as a small- 
scale magneto-convection pattern. At the SST resolution it ap- 
pears similar to that of field-free convection, with upflows in the 
centers and downflows at the edges. However, the correlation 
between continuum intensity and LOS velocity is much weaker 
for this small-scale velocity pattern than for field-free granula- 
tion. This may in part be due to the small horizontal scales of 
these structures combined with the height of formation of the 
630.25 nm line, more easily leading to decorrelations between 
these measured properties than for large (field-free) granules. 
We speculate that the main reason for this decorrelation is that 
small-scale convection does not overshoot to the same heights as 
field- free convection and leaves only very weak traces at the line 
forming and perhaps even continuum forming layers. This would 
imply that the measured RMS velocity is not a good indicator 
of the convective flux for these structures and that possibly there 
are contributions to the velocity signatures from dynamics above 
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the photosphere that is not of convective origin. In this context 
it is important to emphasize that despite the large RMS veloci- 
ties measured for field-free granulation at the formation height 
of the 630.25 nm line, the convective flux in the quiet sun is 
negligibly small 100 km above the photosphere (Nordlund et al. 
2009, Fig. 21). The RMS velocity measured with the 630.25 nm 
line is therefore only a tell-tale of field-free convection actually 
peaking below the visible photosphere and there is no a priori 
reason for expecting the same signature for small-scale magneto- 
convection. 

Given the uncertainty of an appropriate threshold to define 
the mask outlying this small-scale velocity pattern and the likely 
influence of stray light on our estimates of ^los. we conclude 
that there is a transition to small-scale magneto -convection when 
the field covers a sufficiently large area and reaches an average 
strength of 600-800 G. 

3D magneto-convection simulations with similar average 
field strengths over a large FOV (6x6 Mm, corresponding to 
1 15x115 of our pixels) were discussed by Vogler (2004, 2005), 
while Stein et al. (2003) presented 12x12 Mm simulations made 
with an initial average field strength of 400 G. Vogler (2004, 
Fig. 10) discussed snapshots from simulations with average field 
strengths of 400 G and 800 G. By comparison, the largest aver- 
age field strength within a 6x6 Mm box of our data is 600 G. 
The 400 G simulation of Vogler (2004, 2005) shows granules 
with somewhat reduced size but otherwise looking similar to 
field- free granules. The intergranular lanes are wider than for 
field-free convection, contain nearly all the magnetic field and 
show flows of strongly reduced magnitude. The 800 G simu- 
lations show a much more small-scale intensity pattern and the 
velocity map shows small-scale strong upflows superimposed on 
a background with small velocities. Strong downflows are found 
only in thin sheets at the boundaries of the strong magnetic field, 
also constituting the boundaries of small bright granules. The 
map of the vertical magnetic field made from the simulations 
appears diff'erent from our ^los map. The simulated map shows 
very small nearly field-free granules, corresponding to the loca- 
tions of upflows, superimposed on a background of very strong 
(nearly 2 kG) field. These diff'erences may be primarily related 
to the inadequate spatial resolution of the SST data and/or stray- 
light but also due to the diff'erences in height for the synthetic 
map (made at a continuum optical depth T500 = 1) and the 
map from the Helix inversions (representative of approximately 
T500 ~ 0.01-0.1). We are not aware of any publications present- 
ing similar synthetic maps at heights relevant for the formation 
of the 630.25 nm line from simulations with strong average field 
strengths, needed for further interpretation of our data. 

Vogler et al. (2005) calculated the averaged bolometric disk 
center intensity for 6x6 Mm simulation "boxes" with diff'erent 
average field strengths. He found that average field strengths up 
to about 250 G leads to a small enhancement of the bolometric 
intensity by approximately 1%, but for stronger average fields 
the intensity is reduced. At an average field strength of 400 G, 
the reduction in intensity is only 2% but for 800 G average field, 
the disk center intensity is reduced by 10%. The average inten- 
sity within our mask, including the pores, is 98.6% of that for 
the field-free region. Excluding the pores, this average intensity 
is 99.6%. Our magnetic region thus appears much brighter than 
expected from simulations. This may be a coincidental disagree- 
ment or indicate a problem of relaxation of the simulation model, 
its deeper structure or the lower boundary. The fourth panel in 
the top row of Fig. 2 shows a complicated correlation between 
the measured continuum intensity and the strength of the LOS 
magnetic field. For I^losI < 200 G, the intensity decreases with 




Fig. 3: Locations of magnetic bright downflows (coded as white) 
superimposed on a binary map outlining the region with I^losI > 
200 G, shown gray. Tick marks are in units of arcsec. 




Fig. 4: Map of the Stokes V amplitude asymmetry, scaled from 
-0.1 to +0.5 within a mask defined I^losI > 75 G. The FOV 
shown corresponds to the ROI marked in Fig 1 . Tick marks are 
in units of arcsec. 
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increasing field strength. For I^losI larger than 200 G but smaller 
than 600 G, the intensity increases with I^losI- For stronger 
fields the intensity again decreases with I^losI- Vogler et al. 
(2005) investigated similar relations between the disk center in- 
tensity and the vertical component of the magnetic field, using 
simulations with average field strengths of 0, 50, 200, 400, and 
800 G. He found that the intensity increases monotonously with 
field strength for all simulations, i.e., the strongest fields always 
appear to be the brightest. This is explained as an opacity eff'ect, 
shifting the r = 1 level downward into hotter regions, where the 
field strength is higher. This opacity eff'ect is large when the av- 
erage field strength is low and amounts to only a few percent 
when the average field is strong. A similar variation of the field 
strength with temperature (at T500 = 0.1) was found by Shelyag 
et al. (2007) from a 200 G simulation. In our observations, re- 
gions with I^LOsI stronger than 900 G correspond mostly to mi- 
cropores, the formation of which is missing in the simulations of 
Vogler and Shelyag et al. Such pores did however form sponta- 
neously in the 400 G simulations of Stein et al. (2003) and were 
also found in the 200 G simulations of Vogler et al. (2005). 



3.2.1. Magnetic downflows 

The lower three panels in Fig. 2 show the correlation between 
the LOS velocity vlos and I^losI for dark (/c/(/c) < 0.9), aver- 
age (0.9 < /c/(/c) < 1.05 and bright IJ{h) > 1.05) features. 
The dark structures mostly correspond to weak (I^losI < 100 G) 
or strong I^losI > 800 G fields. At locations of very weak 
fields, the dark, average and bright structures show the expected 
signatures of convection: dark downflows and bright upflows. 
However, as I^losI increases, the average LOS velocity increases 
such that all structures show a strong tendency for downflows 
when \Bios\ is in the range 100-400 G. The most striking plot 
is that for the bright structures (third panel in the lower row of 
Fig. 2). At low field strengths, bright structures show average 
upflows of 1 kms"^ but 300-500 G bright structures show av- 
erage downflows of 300 ms"^. Figure 3 shows the locations of 
all downflows (coded as white) stronger than 300 ms"\ hav- 
ing continuum intensities stronger than 1.05 and I^losI > 50 G. 
Such bright magnetic downflows are over 30 times more fre- 
quent per unit area than for I^losI < 50 G. Also shown as gray is 
the region defined by the I^losI > 200 G mask. It is evident that 
virtually all bright magnetic downflows occur near the bound- 
ary to the strong field. At these locations "leakage" of polar- 
ized straylight from the limited spatial resolution and straylight 
is likely to have a significant influence, such that the actual field 
strengths for these features may be significantly lower than mea- 
sured. Our result is qualitatively consistent with the simulations 



of Vogler (2004), showing strong downflows only in thin sheets 
at the interface between field-free granular structures protruding 
into the stronger field. The brightness of these downflows is ex- 
plained by the transparency of the magnetic gas, allowing us to 
see deeper into hotter layers of the downflowing (nearly field- 
free) gas (Keller et al. 2004; Carlsson et al. 2005). The fourth 
plot in the lower row of Fig. 2 shows the relation between vlos 
and I^LOsI based on all data points. This is similar to Figs. 5 and 
6 of Morinaga et al. (2008), also showing the tendency for more 
downward flows as I^losI increases from zero to about 400 G. 
However, the amplitude of the vlos variation is much larger for 
the SST data than for the Hinode data. Apart from the higher 
spatial resolution of the SST data, a difl'erence between the anal- 
ysis of Morinaga et al. and ours is that their LOS velocities were 
obtained by fitting the line core of the 630.15 nm line, thus cor- 
responding to velocities at a higher layer in the atmosphere than 
for our data. Also, Morinaga et al. calibrated their LOS velocities 
to the average value for the surrounding quiet sun (showing con- 
vective blueshift), whereas we use the pores within the FOV as 
reference. Compensating for the diff'erences in zero point, we ob- 
tain qualitatively similar results: an increased LOS velocity for 
field strengths in the range 0-300 G. We also note the presence 
of a weak population with strong downflows for bright structures 
with I^LOsI in the range 200-600 G. 

Finally, we note (third panel in upper row of Fig. 2) that 
both dark and bright features with the LOS magnetic field 
stronger than 600 G on the average show upflows, reaching up 
to -500 m s"^ for the strongest fields. 

3.2.2. Stokes V asymmetries 

To exploit observed properties of the Stokes V profiles ignored in 
Milne-Eddington inversions, we discuss also their asymmetries. 
Several such measures of asymmetries are used in the literature: 
the net circular polarization, the area asymmetry and the am- 
plitude asymmetry. The limited wavelength range scanned with 
CRISP for these data (c.f.. Fig. 7) prevents accurate estimates of 
the first two of these measures. Here, we discuss only the Stokes 
V amplitude asymmetry Sa, defined as 



Sa = (^b - aY)/(a\^ + a^) 



(1) 



where a^ and ar are the absolute values of the maximum and 
minimum values of the Stokes V profile in the blue and red 
wings respectively. These were estimated by fitting second order 
polynomials to V at the three wavelengths closest to the extrema 
and from the fitted coefficients determine max/min values. For 
a small percentage of pixels, the extrema are at the first or last 
wavelength scanned by CRISP; for these we simply used the 



G. Narayan and G.B. Scharmer: Small-scale convection signatures associated with strong plage solar magnetic field 



maximum or minimum values found to estimate da. Since we 
use only three wavelengths for the fits to each of the extrema, 
the measured amplitude asymmetry is noisy for field strengths 
below about 100 G. In Fig. 4 we show da scaled from -0.1 
to +0.5 with a mask to eliminate pixels where the strength of 
the LOS magnetic field is less than 75 G. Above this threshold, 
we find no multi-lobed V profiles among the hundreds of pro- 
files inspected. For nearly all pixels shown, the Stokes V ampli- 
tude asymmetry is positive. Strongly enhanced 6a forms irreg- 
ular connected "networks", shown bright in the figure. Blinking 
the da map with the continuum image, we note that strong Stokes 
V amplitude asymmetry is found almost exclusively at the bound- 
ary between regions with small-scale abnormal granulation and 
much larger "normal-looking'' granules. The left-most plot in 
Fig. 5 shows the correlation between I^losI and 6a and includes 
all pixels in the FOV, the red curve corresponds to binned av- 
erages. Quite clearly, da peaks for the weaker fields and Stokes 
amplitude asymmetries are very small for the strongest fields. 
The second plot shows averaged values for dark (Ic/ih) < 0-9, 
upper curve), and bright structures (Ic/ih) > 1-05, lower curve). 
This demonstrates that darker structures have somewhat more 
asymmetric V profiles than brighter structures at small I^losI- 
The third plot shows that downflows on the average have much 
more asymmetric V profiles than upflows when I^losI is in the 
range 200-800 G, although the scatter is large. The fourth plot 
shows only a weak correlation for I^losI stronger than 800 G, 
but with a clear tendency for some upflows to have negative V 
asymmetries. We note that the scatter in Sa is much smaller for 
large \Bios I than for weak fields, this is not an eff'ect of noise. 

The first plot in Fig. 5 corresponds to Fig. 15 of Shelyag 
et al. (2007), showing the variation of the 630.25 nm Stokes V 
amplitude asymmetry with field strength at r = 0.1 from 250 G 
simulation data. This plot looks strikingly similar to ours, except 
that the simulated data shows an amplitude asymmetry that is 
50% larger for small field strengths than for our data. It appears 
plausible that this diff'erence is related to the lower spatial and/or 
spectral resolution of the CRISP/SST data. 

We note that the Stokes V asymmetry peaks at roughly the 
same low field strengths for which dark as well as bright gas 
on the average shows downflows c.f.. Fig. 2. This strongly sup- 
ports the explanation of the V asymmetries by Shelyag et al. 
(2007) as due to the peripheral parts of magnetic flux concentra- 
tions expanding with height and being surrounded with (strong) 
downflows. At these locations, the LOS cuts through the mag- 
netopause such that the magnetic field strength decreases with 
depth while the downflow velocity increases in magnitude with 
depth, consistent with Sa being positive. Interpretations along 
these lines were first proposed by Sanchez Almeida et al. (1988) 
and supported by, e.g., Grossmann-Doerth et al. (1988); Solanki 
(1989); Martinez Fillet et al. (1997). Our present high-resolution 
data confirms this nicely by locating the peaking of Sa precisely 
at the boundary between abnormal and normal-looking granules, 
where we expect magnetic canopies to expand over the more 
field- free gas. The similar correlations for weak plage simulation 
data and the present observed data from strong plage suggests 
that the dynamics close to the boundary of flux concentrations is 
similar for small and more extended flux concentrations. 

The fourth plot in Fig. 5 shows the correlation between V 
asymmetry and LOS velocity for field strengths stronger than 
800 G, corresponding to the interior of the flux concentrations. 
Here, the gas shows mostly upflows (see lower four plots of 
Fig. 2) and we expect the field to expand with height and its 
strength thus to decrease with height. In this region Sa remains 
mostly positive, as expected with upflow velocities decreasing 



in magnitude with height. The correlation between vlos and Sa 
is weak, but in the sense of being stronger for downflows than 
for upflows, for pixels with I^losI stronger than 800 G. There is 
also a small fraction of pixels with negative V asymmetry where 
upflows are stronger than 500 ms"^ Whether this fits within a 
magneto-convection scenario with strong field remains to be in- 
vestigated. 



3.3. Small-scale velocity features 

Figures 6 and 9 shows the ROI chosen for identification and dis- 
cussion of specific velocity features. The two upper panels show 
the continuum intensity and line minimum intensity /min in the 
630.2 nm line. Comparing to Fig. 1 and the ^los map in the 
upper-right panel of Fig. 6 indicates that the line minimum inten- 
sity is a good proxy for detecting strong-field magnetic features. 
This is in part because of line weakening due to the larger tem- 
perature at equal optical depth in the magnetic flux concentra- 
tions (Chapman & Sheeley 1968; Shelyag et al. 2007). The vlos 
map, made from data returned by Helix is shown in the lower- 
right panel of Fig. 6. The white/black/red contours in these pan- 
els correspond to I^losI = 200 G and outline the same regions 
as high-lighted in the Vlos map of Fig. 1. 

The line of sight velocity (vlos) map in Fig. 6 shows the 
signatures of field- free and magneto-convection discussed ear- 
lier. There are a few areas with strong downflows reaching up to 
3.5 km s"^ ; the strongest upflows are approximately -1.3 km s"^ . 

The bright features in the line minimum images appear suf- 
ficiently similar to those in images recorded in the G-band and 
the Ca II H line to allow us to describe properties of the observed 
features following the nomenclature of Berger et al. (2004) and 
Rouppe van der Voort et al. (2005). In the following, we discuss 
the results of the inversions for a few selected features of each 
type, indicated with circles in Fig. 6. For each of these selected 
features, we show the observed and fitted Stokes / and V pro- 
files at the center of the corresponding circle (Figs. 7 and 10) 
and the inferred variation of vlos and ^los along a horizontal 
cut through the center of the same circle (Figs. 8 and 11). 

In Figs. 7 and 10 we show observed and fitted Stokes / and V 
profiles at center of the features discussed in the following sub- 
sections. Fits for which the Stokes V signal is strong appear good 
and justify the use of Milne-Eddington inversions neglecting 
LOS gradients of the magnetic field and LOS velocity. However, 
for the profiles with the weakest V signal shown in Fig. 10, the 
fits are quite poor. 



3.3.1. Pores 

The strongest flux concentrations appear as dark pores in the I^ 
image but are only marginally darker than their surroundings in 
the /min map. In the following, we refer to tiny pores which ap- 
pear as distinctly dark in the continuum intensity (I^) maps as 
micropores and those that have only hints of darkening in the 
continuum intensity, but still show strong flux concentrations in 
the ^Los map, as protopores. With this definition, the feature 
labeled M in Fig. 6 is a micropore and the one labeled P is a 
protopore. Most of the pores show at least one strong down- 
flow channel adjacent to their boundaries. This is in agreement 
with earlier observational studies, e.g. by Leka & Steiner (2001), 
Sankarasubramanian & Rimmele (2003), Hirzberger (2003) and 
Stangl & Hirzberger (2005). Thin downflows lanes around pores 
were also seen in the simulations of Cameron et al. (2007). 
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Fig. 6: The ROI shown in Fig. 1, used to identify isolated bright points (B0-B3), ribbons (R0-R3), flowers (F0-F3), strings (S0-S3) 
and other features (00-03). The upper row shows the continuum image and the /min map, the lower row ^los and vlos maps 
obtained from inversions. The white/black/red contours correspond to I^losI = 200 G. Tick marks are in units of arcsec. The color 
bar indicates the signed LOS magnetic field in Gauss, the grey scale bar the LOS velocity in ms"^ 



3.3.2. Isolated bright points 

Circles labeled BO, Bl etc in Fig. 6 contain examples of isolated 
bright points. Figure 8 shows the plots of vlqs and ^los across 



the various features. The top row of Fig. 8 show plots across the 
bright points. For the feature marked BO there is a strong down- 
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Fig. 7: Observed (solid) and fitted (dashed) Stokes / and V profiles for isolated bright points (B0-B3), ribbons (R0-R3), flowers 
(F0-F3), strings (S0-S3) and other features (00-03) indicated in Fig. 6. Both Stokes / and V are normalized to the continuum 
intensity, the wavelength of which is outside the range plotted. Wavelengths along the x-axis are in units of mA. 



flow of about 960 ms"^ close to where the field is strongest. 
In B 1 we can also see that there is a strong downflow of about 
2.3 kms"^ close to where the LOS field peaks. In B2 we see a 
similar though weaker downflow of 840 m s"^ . We note that the 
downflow occurs adjacent to but not precisely where the LOS 
field peaks. This is consistent with previous studies of G-band 
bright points (Berger et al. 2004; Rimmele 2004). B3 seems to 
be an exception with no strong downflow but instead an upflow 
of about -430 ms"^ close to the location of the strongest LOS 
field. 



3.3.3. Ribbons and flowers 

The micropores are surrounded by several small scale magnetic 
features that form almost circular interfaces between the microp- 
ores and (the disturbed) surrounding granules. These small scale 
features are clearly visible in the line minimum images as bright 
structures having ribbon- and flower-like shapes. Similar fea- 
tures were first observed and described by Berger et al. (2004) 
and further studied by Rouppe van der Voort et al. (2005). 
Ribbons and flowers are similar structures, the difference being 
that flowers are somewhat circular and ribbons are rather elon- 



gated. When ribbons take a circular shape, they are here referred 
to as flowers. The features marked F0-F3 etc are examples of 
flowers and those marked R0-R3 are examples of ribbons. The 
/min maps shows that the ribbons and flowers have somewhat 
dark cores and brighter edges. The ^los map confirms that the 
flowers and ribbons are magnetic. An aspect worth noting is that 
in our observations flowers are found only in the close vicin- 
ity of micropores and protopores. Ribbons are found mostly in 
the vicinity of micropores and protopores and in nearby network 
plage regions. 

The second row of Fig. 8 shows plots across ribbons. RO 
shows a strong upflow of about -840 ms"^ very close to the 
peaking of the field. This is also seen in Rl where the upflow 
peaks at -450 ms"^. The region of strong ^los seems slightly 
more extended in Rl. Although we see similar variations in R2, 
the upflow is much weaker reaching only -290 m s"^ . The region 
of strong ^Los seems to be slightly more extended for R3 than 
for R2. We see a primary upflow of -420 m s"^ within the feature 
and a secondary upflow of -680 m s"\ in contrast to the single 
upflows found for other ribbons. 

The third row of Fig. 8 shows plots across flowers. FO shows 
a strong upflow of about -760 ms"^ close to the peaking of 
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Fig. 8: Horizontal cuts of LOS velocity (dashed) and strength of LOS magnetic field I^losI (solid) for isolated bright points (BO- 
BS), ribbons (R0-R3), flowers (F0-F3), strings (S0-S3) and other features (00-03) indicated in Fig. 6. Distances along the x-axis 
are in units of arcsec. 



the field. Fl is similar but with an upflow of -580 ms"^ F2 
shows similar variations but with a smaller upflow of around 
-400 ms"^ F3 although shows a upflow close to -400 ms"^ 
seems to have difl'erent properties compared to the other flower 
features. This could be because it is at the boundary of a fully 
developed micropore and hence we do not see the peaking of the 
field at the flower as clearly. The extended nature of flowers is 
clear from these plots. Flowers show only a small variation ^los 
but strong variation in vlos • Another common feature is that at 
the boundary of the flowers there is often a transition from a 
downflow to an upflow. 



3.3.4. Strings 

The network also consists of other bright features such as thread- 
like features classified as strings by Rouppe van der Voort et al. 
(2005). The features labeled SO, SI etc are example of strings. 

The fourth row of Fig. 8 shows plots across strings. In SO we 
see a downflow of about 300 m s"^ close to the peak field. In SI 
we see a downflow of about 120 ms"^ In S2 we see a strong 
downflow of 1.4 kms"^ S3 is difl'erent from the other strings 
studied and shows an upflow of about -640 ms"^ close to the 
position of the peak LOS field. There is also a secondary upflow 
of about -1.0 kms"^ nearby. This secondary upflow appears to 
be related to a tiny bright structure visible in the continuum map. 
It is not clear if this structure is a tiny fragment from the neigh- 
boring granule, ^los shows that it is a weak magnetic structure. 
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Fig. 9: The ROI shown in Fig. 1, used to identify various upflow features (marked A-T). The upper row shows the continuum image 
and the /min map, the lower row ^los and vlos maps obtained from inversions. The black/red contours correspond to I^losI = 200 G. 
Tick marks are in units of arcsec. The color bars indicates the signed LOS magnetic field in Gauss (left) and LOS velocity in 
m s"^ (right). Except for the color and grey scales used, this figure is identical to Fig. 6. 



3.3.5. Other small-scale magnetic features 

The fifth row of Fig. 8 shows plots of a few miscellaneous small- 
scale magnetic features which we refer to as "other" features. 



These are bright features which do not fit in any of the above 
categories. Some examples are shown in circles labeled as 01, 
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02 etc. OO looks like a fragmented string and shows a strong up- 
flow of about -1.05 kms"^ 01 shows a variable downflow pat- 
tern ranging between 100 m s"^ and 600 m s"^ 02 a similar flow 
pattern in the range 200-800 ms"^ 03 is an extended structure 
with an upflow varying between -300 m s"^ and -450 m s"^ 



3.3.6. Small-scale upflow features 

From the vlos map in Fig. 9 it is clear that there are sev- 
eral small-scale upflows often surrounded by ring-shaped down- 
flows. In order to study these we randomly identified and se- 
lected strong upflows using the vlos map. These upflow features 
are marked with circles labeled A-T. 

Features A-D, N and R are all small scale upflows occurring 
in regions where I^losI is in the range 10-230 G. The upflows, 
based primarily on the Stokes / profiles, for all these cases are 
fairly strong, ranging from -900 ms"^ to -1.4 kms"^ The sur- 
rounding downflows range between 550 ms"^ and 1 kms"^ In 
cases A, B, D, N and R we see a weakening in the LOS mag- 
netic field at the location of the upflow. In C we see the ^los 
drop off, the reason being its location at the interface of a strong 
field and a field-free region. From the corresponding I^ image it 
appears that A, B, D, and R are probably granular fragments. N 
and C seem to be located in intergranular lanes but it is difficult 
to label these features from the I^ image. In the /min map, N and 
C appear completely dark and they are obviously different from 
the brighter features A, B and D in this respect. We note, that for 
some of these features (especially B, but also N and R), the ob- 
served Stokes V profiles appear strongly redshifted with respect 
to the fitted profiles. For feature B, this redshift is approximately 
2 km s"^ such that the Stokes V profile actually indicates a down- 
flow of about 500 ms"^ The very weak Stokes V profile mea- 
sured, peaking at less than 0.7%, indicates the possibility that 
we actually measure polarized straylight from the neighboring 
magnetic feature. 

Features F-H, J-L and S represent small-scale upflows in 
regions of where I^los I is in the range 60 G to 1.1 kG. For these 
features we see upflows ranging from -500 m s"^ to - 1 .0 km s"^ 
The surrounding downflows range from 200 m s"^ to 1.1 km s"^ 
In all of these cases, with the exception of H, we see that the 
peaking of Vlos coincides with peaking of ^los- F appears to be 
part of the ribbon RO of Fig. 6 and G a part of flower FO. J is part 
of a flower-like feature and L is a part of a ribbon-like feature. 
K happens to be a part of OO feature of Fig. 6. H seems to be a 
fragment of a granule which is adjacent to a ribbon. 

Features marked E, I, M, O and P correspond to upflows in 
regions of I^losI ranging from 5 G to 650 G. The upflows range 
from -530 ms"^ to -1.45 kms"^ The surrounding downflows 
are in the range 550 m s"^ to 700 m s"^ E, I, O, and M appear to 
be fragments of granules. P appears to be an isolated bright point 
with an upflow at its center and a strong downflow of 1.4 km s"^ 
at its boundary. The field peaks at the center where the upflow is 
located. 

Finally we have feature like Q and T which are located in 
very weak field regions with ^los ranging between and 40 G. 
The upflows range between -550 ms"^ and -900 ms"^ The 
surrounding downflows are in the range 550 m s"^ to 1.0 km s"^ 
These are most likely granular fragments. 

From this we have gained additional support for flowers and 
ribbons being related to small-scale upflows. Fig. 9 shows the 
observed and fitted stokes profiles for A-T. 



3.3.7. Discussion 

From these plots we conclude that the various features we have 
described have specific signatures in ^los and vlos- Ribbons 
and flowers seem to possess strong upflows while isolated bright 
points are associated with strong downflows. Strings, on the 
other hand seem to show downflows of intermediate strength. 
We have shown only a few examples for each type here but have 
examined also several other features in the data. Although some 
exceptions can be found, the differences in flow patterns de- 
scribed for these features appear systematic. We have also found 
various small-scale distinct features which do not fit into any of 
the categories defined previously. 

There have been a few recent studies of magnetic upflows 
and downflows in magnetic elements at 0'.'l-0'.'2 resolution by 
Berger et al. (2004) and Rouppe van der Voort et al. (2005). 
Upflows were found within a few ribbons studied by Rouppe 
van der Voort et al. (2005). Langangen et al. (2007) made a spec- 
troscopic study of magnetic features such as flowers and found 
weak upflows of about 150 m s"^ in their centers, similar to what 
we find here. Rimmele (2004) studied magnetic fine structure 
in a region close to disk center based on magnetograms in the 
blue wing of the Fe i 630.25 nm and LOS velocities estimated 
from filtergrams in the blue and red wings of the weak C i line at 
538 nm and the Fei line at 557.6 nm. Without making the same 
detailed distinction between different types of magnetic structure 
as made here, he found downflows in the range a few hundred 
m s"^ up to 1 km s"\ located at the edges of small flux concentra- 
tions and chains of bright points. He also found that the down- 
flows appear significantly narrower and faster at deeper layers 
and very weak flows within the flux concentrations at heights 
corresponding to the 630.25 and 557.6 nm lines. The average 
velocities within the bright structures were found to be close to 
zero at higher layers and a few hundred m s"^ in the deeper lay- 
ers. 

The ubiquitous presence of downflows near the boundaries 
of strong fields found in the present study appear consistent with 
the simulations of Vogler et al. (2005), showing stronger such 
downflows. In other simulations made at different average field 
strengths (Vogler 2004), velocity maps were made for heights 
corresponding to T500 = 1 • The morphology of our velocity maps 
look like a mixture of the simulated velocity maps obtained with 
average field strengths of 5 = 400 G and 5 = 800 G. 



4. Conclusion 

In this paper, we analyze spectropolarimetric S ST/CRISP data, 
obtained at a spatial resolution close to 0^.^15, from a unipo- 
lar ephemeral region with small pores close to sun center. We 
use MOMFBD methods (Lofdahl 2002; van Noort et al. 2005) 
to restore and align the polarimetric images recorded at differ- 
ent wavelengths and the Milne-Eddington inversion code Helix 
(Lagg et al. 2004) to retrieve LOS velocities and the LOS com- 
ponent of the magnetic field. These inversions were made with 
the magnetic filling factor set to unity, such that LOS magnetic 
fields obtained and discussed correspond to measured flux den- 
sities. 

Using this data, we investigate signatures of magneto- 
convection associated with strong average magnetic fields. We 
use a 200 G threshold for the LOS magnetic field to define 
such regions of strong field but emphasize that the average 
field strength within them is approximately 600-800 G, the up- 
per limit taking into account a likely underestimate of the field 
strength due to our neglect of spatial straylight. This data corre- 
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sponds to unipolar strong field over extended contiguous regions 
of typically T' width or more across the smallest dimension. 
The field strength within a larger 6x6 Mm box, used in MHD 
simulations relevant to the present work, reaches up to at least 
600 G. We discuss statistical relations between measured inten- 
sities, LOS velocities and LOS magnetic fields. In addition, we 
identify specific structures such as pores, bright points, strings, 
ribbons and flowers (Berger et al. 2004) and discuss their LOS 
velocity and magnetic field signatures. We compare our results 
to earlier analysis of similar data and to 3D MHD simulations. 

Within the strong flux concentrations, we find a small-scale 
granular velocity pattern resembling that of field-free granules 
but with scales about 4 times smaller. Based on a zero-point for 
LOS velocities determined by assuming the pores to be at rest, 
this small-scale velocity field averages to nearly zero (70 m s"^) 
and mostly corresponds to upflows at the centers of magnetic 
"granules", surrounded by downflows. These LOS velocities 
correlate, though weakly, with the continuum intensity such 
that upflows are on the average brighter than downflows. The 
weak correlation between continuum intensity and LOS velocity 
found within the strong flux concentrations, compared to what 
is the case for field-free convection, could be taken as evidence 



against the interpretation of a convective origin for this velocity 
field. However, a convective flux below the visible surface is ob- 
viously needed to explain the high radiative energy flux in these 
extended (approximately 2'' across the smallest dimension) mag- 
netic structures. This suggests to us that the main reason for the 
weak correlation rather is that the nature of small-scale magneto- 
convection is distinctly difl'erent, in particular as regards leav- 
ing "tell-tales" of convection in the line-forming layers above 
the photosphere, from that of large-scale field-free convection. 
Given the uncertainty of an appropriate threshold to define the 
mask outlying this small-scale velocity pattern and the likely in- 
fluence of stray light on our estimates of ^los^ we therefore con- 
clude that there is a transition to small-scale magneto-convection 
when the field covers a sufficiently large area and reaches an av- 
erage strength of 600-800 G. The measured RMS velocity of this 
magneto-convection is 400-500 ms"^ compared to 700 ms"^ 
for the nearly field-free surrounding granulation. The measured 
RMS velocity decreases with increasing field strength for dark 
as well as bright structures. For I^losI > 400 G, the signed 
average LOS velocity decreases systematically with increasing 
I^LOsI such that dark and bright structures on the average show 
upflows when I^losI > 700 G. The RMS continuum intensity 
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Fig. 11: Horizontal cuts of LOS velocity (dashed) and strength of LOS magnetic field I^losI (solid) for upflow features (A-T) 
indicated in Fig. 9. Distances along the x-axis are in units of arcsec. 



variations within the mask are nearly as strong inside the mask 
as outside the mask. 

The boundaries of the magnetic regions, defined in terms of 
a LOS field of about 200 G, mark the transition between normal- 
looking granules and small-scale (abnormal) granulation. Here, 
we find predominantly downflows. In particular, there is a pop- 
ulation of bright downflows that are over 30 times more fre- 
quent than in a field- free area of the same size. At this mag- 
netic boundary, the Stokes V profiles are asymmetric with the 
blue peak stronger than the red peak, in agreement with syn- 
thetic y profiles calculated from simulations with 250 G aver- 
age field strength (Shelyag et al. 2007). We identify this popula- 
tion with strong downflows occurring at the interfaces between 
strong magnetic field and field-free granules, seen in simulations 
also with strong (400-800 G) average field Vogler (2004), re- 
producing the small-scale magnetic granulation pattern seen in 



our data. The similarity of the correlation plots shown in Fig. 2, 
with contributions mostly from extended regions of strong field, 
and those of Rimmele (2004), made from selected bright points 
and other small-scale flux concentrations, as well as similar plots 
made from 250 G simulations (Shelyag et al. 2007), suggests 
that the dynamics close to their boundaries is similar for small 
as well as large flux concentrations. 

The weak field returned by Helix close to the boundaries of 
the flux concentrations can be interpreted as partly due to mag- 
netic canopies expanding over the more field- free adjacent pho- 
tosphere with normal-looking granules. The gradually weaker 
field away from the boundaries obtained from the inversions is 
then not only from polarized straylight but also an eff'ect of the 
magnetopause gradually shifting to smaller optical depths away 
from the flux concentration, as discussed by e.g., Martinez Pillet 
etal. (1997). 
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A detailed comparison with our results and those from simu- 
lations as regards the interior parts of flux concentrations, where 
the field is strong over extended regions, has not been possible. 
This requires the quantities returned by Helix to be compared to 
simulation data at an optical depth T500 = 0.01-0. 1, but presently 
only simulations with strong magnetic field have been analyzed 
at a height corresponding to T500 = 1 • 

We also compare our correlation plots between the LOS ve- 
locity and LOS magnetic field with measurements made at lower 
spatial resolution based on HINODE/SOT spectropolarimetric 
data (Morinaga et al. 2008). We found good qualitative agree- 
ment as to the variation of the average and RMS LOS velocity 
with field strength, but the RMS velocities measured with the 
present data are much higher than measured with Hinode. To 
some extent this diff'erence can be attributed to diff'erences in the 
methods used to measure LOS velocities - Morinaga et al. use 
fits to the line core of the 630.15 nm line, whereas we use the 
less opaque 630.25 nm line and velocities fitted by Helix to the 
entire line profiles, representing deeper layers. It is however ob- 
vious that the higher spatial resolution of the SST is crucial for 
resolving these very small-scale magnetic velocity fields. 

Using the line minimum intensity map, we identify bright 
points, strings, ribbons and flowers (Berger et al. 2004) and dis- 
cuss their LOS velocity and magnetic field signatures. We find 
preferred upflows in ribbons and flowers while isolated bright 
points and strings mostly show downflows. This is in good agree- 
ment with earlier studies, e.g. (Rouppe van der Voort et al. 2005) 
and (Langangen et al. 2007). Most of these specific structures 
identified are located within the 200 G mask used to outline the 
small-scale magneto-convection pattern seen in the LOS veloc- 
ity map; they are thus part of that convection. 

The data presented clearly demonstrate the importance of 
high spatial resolution: the convection features shown in the 
present S ST/CRISP data are barely resolved and would be dif- 
ficult, if not impossible, to resolve with a significantly smaller 
telescope. A somewhat surprising shortcoming is the lack of 
simulation data with strong average magnetic fields (400-800 G) 
and properties extracted from heights relevant to the formation 
of the 630.25 nm iron line. The simulation results discussed by 
Vogler (2004) are presented only for T500 = 1 • Such future com- 
parisons obviously are desirable and will need to take into ac- 
count the limited spatial resolution as well as straylight of the 
telescope used to record the data. We finally note that the Stokes 
Q and U profiles, although mostly weak, show up clearly at some 
locations, allowing additional consistency tests such as investi- 
gations of their asymmetries relative to those of the Stokes V 
profiles (Martinez Fillet et al. 1997, Appendix B). 
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